Multiple Avalanches Across the Metal-Insulator Transition of Vanadium Oxide 
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The metal insulator transition of nano-scaled VO2 devices is drastically different from the smooth 
transport curves generally reported. The temperature driven transition occurs through a series of 
resistance jumps ranging over 2 decades in amplitude, indicating that the transition is caused by 
avalanches. We find a power law distribution of the jump amplitudes, demonstrating an inherent 
property of the VO2 films. We report a surprising relation between jump amplitude and device size. 
A percolation model captures the general transport behavior, but cannot account for the statistical 
behavior. 

PACS numbers: 71.30.+h, 72.80.Ga, 64.60.Ht,64.60.an 
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There are many systems in nature that have a tran- 
sition from one state to another which is driven by an 
external force, where the transition is not continuous, 
but rather through a scries of avalanches. These sys- 
tems are diverse as is the driving force and the measure- 
ment technique used to observe them. Examples include: 
Barkhausen noise in ferromagnets p], 0] , acoustic emis- 



, magnetocaloric effect 
sharp magnetization 



sion in martensitic transition 
in giant magnetocaloric alloys 

steps or sharp resistance steps in manganites @, 0] and 
capillary condensation of He in nanoporous material Q . 
In general an external parameter modifies the free energy 
of the two phases, which provides the driving force for 
the system. The nature of the avalanches provides much 
information about the system at hand. One can learn 
about the types of interaction, the role of fluctuations, 
existence of self organized criticality and the universal 
features of the transitions which transcend the specific 
physical system H 0, EE E I ■ 

Several of these systems are characterized by phase 
separation during the transition. This implies that the 
transition occurs through a series of avalanches trans- 
forming portions of the system from one phase to the 



KJP 

other 12j, ll3j, ll4( . One such system, which has received 
much attention, is Vanadium Oxide (VO2). VO2 un- 
dergoes a first order metal insulator transition (MIT) 
of over 4 orders of magnitude at ~ 340 K. The transi- 
tion is from a high temperature metallic rutile phase to a 
low temperature insulating monoclinic phase, and can be 
driven by temperature, light irradiation or pressure [l5j |. 
Avalanches may be expected in this system for multiple 
reasons; it has a first order phase transition, there is a 
state of phase separation between metallic and insulating 
regions along the transition 3, 17 1; and ultra fast mea- 
surements reveal a phase transition in separated domains 
of the system with a transition time on the order of a few 
picoseconds 18, In this paper we report the first ob- 
servation of multiple avalanches across the temperature 
driven MIT in V0 2 . 



Generally, in order to identify an avalanche, the sen- 
sitivity of measurement has to be greater than the mag- 
nitude change of the relevant parameter, and the mea- 
surement frequency faster than the avalanche frequency 
0- 0, [2il • For resistance measurements in VO2 this im- 
plies that the size of the device has to be comparable 
to the magnitude of the domains involved in a single 
event, and that the rate of change of the driving force 
slow enough to resolve the avalanches. 

By measuring transport properties of VO2 devices, 
with electrode spacing as small as 200 nm and tempera- 
ture sweeping rates smaller than 3K/min, we were able 
to detect that the MIT evolves through a series of dis- 
creet jumps ranging over two decades of resistance in- 
dicating the transition is through a series of avalanches. 
We find a scaling law which shows that the amplitudes of 
the largest jumps are inversely proportional to the device 
length. Interestingly the jump size distribution follows a 
power law which is generally robust for various impor- 
tant parameters, including different depositions, device 
size and temperature sweep rates. The general shape 
of the transport measurement can be understood in the 
framework of a non interacting site percolation model. 
But it cannot explain the statistical nature of the jump 
sizes, implying that the power law distribution is an ex- 
pression of an intrinsic property of the VO2 films. 

Vanadium oxide thin films were prepared by reactive 
RF magnetron sputtering of a vanadium target (1.5" di- 
ameter, > 99.8%, ACI Alloys, Inc.) on an r-cut sapphire 
substrate. The samples were prepared in a high vac- 
uum deposition system with a base pressure of 5 x 10 -8 
torr. A mixture of ultra high purity (UHP) Argon and 
UHP oxygen gasses were used for sputtering. The to- 
tal pressure during deposition was 3 x 10~ 3 torr, and the 
oxygen partial pressure was optimized to 1.5 x 10 -4 torr 
(5% of the total pressure). The substrate temperature 
during deposition was 500°C while the RF magnetron 
power was kept at 300 W. These conditions yielded a de- 
position rate of 0.37 A/s, and a total thickness of 90 nm 
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FIG. 1: Main panel- 8 consecutive R-T cycles of a 0.5x6 \imr 
VO2 device zoomed in on the area marked in the full mea- 
surement shown in b. a- Image of one sample with device 
length varying between 1 to 6 pan and device width of 6 fim. 



FIG. 2: Average value of each of the three largest jumps as a 
function of electrode spacing. The inset plots only the aver- 
age of the largest jump as a function of inverse the electrode 
spacing. Error bars are standard deviations. 



was deposited for all the samples reported here. The 
samples were cooled at a rate of 13°C/min in the same 
Ar/02 flow of the deposition. Films were characterized 
and verified to be single phase VO2 by x-ray diffraction 
(Rigaku RU-200B diffractometer) using Cu Ka radia- 
tion and Energy-dispersive X-ray spectroscopy (Phillips 
XL30 ESEM). Macroscopic measurements of resistance 
vs. temperature (R-T) curves reveal an MIT starting at 
~ 340 K with a change of almost 4 orders of magnitude 
in the resistance, indicative of high purity thin films 21 1. 
The VO2 nano-structures were fabricated by standard e- 
beam lithography and lift-off techniques. On top of the 
VO2 film metallic electrodes were deposited with 200 nm 
to 4 micron separation and 2 um to 15 fj,m width. The 
electrodes were deposited by sputtering of 50 nm Vana- 
dium, which acts as an adhesion layer, and 100 nm of gold 
to assure that the electrode resistance is lower than that 
of the metallic VO-i- A typical sample with 8 devices is 
depicted in of Fig.QJ,. Using standard photolithography, 
the rest of the VO2 was etched in argon plasma, leaving 
only a square of VO2 (clearly seen in Fig.QJ,). In the final 
step the electrodes were photolithographically connected 
to macroscopic pads. The R-T characteristics were ac- 
quired in a home made test-bed by standard 2 probe or 
quasi 4 probe measurements (with no significant differ- 
ence between the two) using a constant current source 
at a measurement rate of 10 Hz. Temperature was con- 
trolled by a Lakeshore 332 temperature controller, with 
sweep rates varying between 0.1 K/min and 3K/min. 
Each measurement was cycled multiple times (10-100) in 
order to obtain sufficient data for statistical analysis of 
the jumps. 

Figure [1] (main panel) shows 8 consecutive cycles of a 
typical R-T measurement across a 1 x 6 urn 2 device fo- 



cusing on the numerous resistance jumps. The full range 
measurement is depicted in Fig. [Tp. The transition tem- 
perature is 338 K and the total change in resistance is 
over three orders of magnitude occurring along ~ 8 K 
with a typical hysteresis of IK, which is similar to the 
macroscopic properties of the sample. The jumps range 
from over 2000 f2 for the largest jumps to below 10 O, 
limited only by the measurement noise. We find that the 
jumps occur between two consecutive measurements even 
for the slowest temperature sweep rates, where the (nom- 
inal) temperature change between adjacent data points is 
smaller than 0.5 mK. This indicates that the time scale 
of each jump is much shorter than our measurement ca- 
pabilities, which is expected due to the fast nature of the 
VO2 transition, of a few picoseconds 18, Tj^. This im- 
plies that the transition between the spatially separated 
but co-existing metallic and insulating phases in VO2 
0, 12 1 is not continuous, but rather occurs through a 
series of avalanches, where each avalanche is expressed 
as a resistance jump in the measurement. 

The main characteristics of the R-T measurement are 
as follow: the first jumps appear close to the onset tem- 
perature of the macroscopic transition. There arc 1-3 
large jumps (ranging up to a few thousands of ohms in 
some cases), which may account for 50% of the resis- 
tance change. The rest of the jumps are smaller, from 
a few hundred ohms and are limited by the resolution 
of our measurement. In some samples, jumps as small 
as 2 f2 were resolvable. Figure [2] shows the three largest 
jumps (each one averaged over at least 10 cycles and 2 
devices) as a function of the device length. Interestingly, 
the largest jump decreases with device length and is a 
linear function of inverse the device length, as demon- 
strated in the inset of Fig.[H We note that the data can 
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FIG. 3: Histogram of the jump amplitudes plotted on a log- 
log scale, acquired from over 100 cycles measured for a sample 
1x6 fim 2 . Solid line (red online) is a linear fit, indicating power 
law dependence, with exponent value a — 2.45. Inset- value 
of a for different channel lengths and different samples. 



be fitted also to an exponential dependence between the 
largest jump and the device length, but not as well. By 
extrapolating the linear relation presented in the inset of 
Fig-Hi to a jump size of zero we find that the R-T would 
look smooth for electrode spacing of over 6 microns. This 
may explain why avalanches were not observed in other 
experiments performed on devices longer than 10 fj,m. 

The interesting nature of the avalanches is illustrated 
by the size distribution of the jumps measured along the 
phase transition. Figure [3] shows a histogram of the am- 
plitudes of the jumps, taken from over 100 cycles, mea- 
sured for a 1 fj.m device at a rate of 1.2 K/min. The 
histogram fits well to a power law, p(A) ~ A~ a , where p 
is a probability function and A is the Amplitude of the 
jump. The value of a from the linear fit is found to be 
2.45 ±0.1. The inset of Fig. [3] depicts exponent values for 
different electrode spacing. Data for the 1 pm device is 
provided from 3 different samples deposited separately, 
but under the same nominal conditions. For these the 
value of a changes between 2.35 and 2.65. In addition we 
find that the power law behavior is identical for heating 
and cooling, and is also not dependant on the tempera- 
ture ramp rate, presented in the inset of Fig.[4ji. To a 
first approximation, the value of the exponent is very ro- 
bust in all these samples, and is not strongly affected by 
any of the parameters tested. Thus, this measured ex- 
ponent is an inherent and characteristic property of the 
VO2 thin films. 

We find evidence for an athermal transition, i.e. ther- 
mal fluctuations do not play an important role in the 
MIT. For one, as pointed out earlier, does not vary much 
as a function of the temperature ramp rate fiol ]. In ad- 



dition we monitored the resistance at constant temper- 
atures along the transition. In most cases, like the one 
presented in Fig. 2k,, there are no jumps when the tem- 
perature is constant. This figure shows that while the 
temperature is ramped at a rate of 0.2 K / minute for the 
first 5 minutes (top curve and right axis) there are a num- 
ber of avalanches and large changes to the resistances 
(bottom line and left axis). But once the temperature 
is stable, during the next 4 hours of measurement there 
are no avalanches. There were a few cases where we did 
observe a single event even though the temperature was 
stabilized. 

The general characteristics we find in our samples can 
be interpreted in context of the percolation nature of the 
measurement system. The observed few large jumps and 
many smaller jumps with a wide range can be understood 
in the framework of a percolation network subject to two 
constraints. First, the resistivity ratio between the insu- 
lating and metallic phases is large but finite (in our case 
Pi I Pm ~ 1000) and second, the area measured is finite, 
so the total number of domains which transit through 
an avalanche is not very large (on the order of pi/pm)- 
In order to model the effects of percolation in the ex- 
perimental system we preformed numerical simulations 
of site percolation on a square lattice. This is reason- 
able, since the macroscopic VO2 MIT can be modeled 
by an effective medium approximation 22( . In our simu- 
lation each site changes randomly from a finite metallic 
resistance to a finite insulating one and the resistances of 
all the sites are identical. Figure [4)d portrays a simula- 
tion with conditions similar to the experimental values. 
The thick (blue online) line is the simulation, while the 
thinner one is an experimental R-T. It is evident that this 
simple toy model captures the main features of the exper- 
imental data. We can identify with certainty the largest 
jump in the experiment as the percolation threshold. If 
we look at the amplitude of the percolation jump for dif- 
ferent device geometries, averaged over 200 realizations, 
we find that the percolation jump does scale with inverse 
the device length in a similar fashion to the experiment 
(shown in the inset of Fig.[4)D). But, we could not re- 
construct the effect of the jump size approaching zero for 
longer devices. 

In the attempt to understand the statistics of the 
avalanches, it is important to note that the amplitude of 
a resistance jump cannot be seen as a direct measure of 
the VO2 volume involved in an avalanche event. Rather, 
this magnitude is superimposed on the percolation na- 
ture of the measurement. This is evident from the wide 
distribution of jumps seen in the simulations (Fig. 3)3) 
even though we assumed that all the sites have identi- 
cal properties. We investigated the contribution of only 
percolation to the power law distribution of avalanches 
using the model mentioned above. In order to do that 
we attempted to calculate values of for device geometries 
similar to the experimental one and in the same fashion 
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FIG. 4: a- Resistance (bottom curve, left axis) and tem- 
perature (top curve, right axis) as a function of time. No 
avalanches are observed for constant temperatures. Inset- a 
for different temperature sweep rates, b- Simulation of R-T 
(thick line, blue online) along with experimental data (thin 
line, black) from a 1x6 /jm 2 device. The simulation is of size 
20x50 (length x width) on a square lattice and pi / pm = 1000. 

as we did for the experimental data. Only for part of the 
data range we could fit the numerical simulations to a 
power law, since as the jumps become smaller the num- 
ber of jumps becomes constant and the slope levels. In 
this range, the value of a varied considerably, between 
1.5 and 3.3, with simulation size and for different his- 
togram binning choices. This means that the distribu- 
tion function of the numerical data is not a power law, 
and the experimental data cannot be explained solely by 
the geometrical nature of percolation. In order to ac- 
count for the distribution of avalanches one has to take 
under consideration intricate conditions of the VO2 sys- 
tem, interactions, inhomogeneities in the sample or in 
the driving force. This work is beyond the scope of the 
current paper. 

In summary, by choosing a properly nanostructured 
film we were able to observe discrete transitions in the 
Metal Insulator transition of VO2 films. These results 
imply that the metal insulator transition occurs through 
a series of avalanches. The largest of these discontinu- 



ous transitions has an inverse dependence with sample 
size. Interestingly, the distribution of resistance jump 
amplitudes follows a reasonably robust (i.c independent 
of experimental parameters) power law with an exponent 
close to 2.4. This value is similar to the ones reported for 
acoustic emission measurements of avalanches in marten- 
sitic transitions 0, E3] j however in our case the measured 
amplitude is also a consequence of the percolating nature 
of the system. Importantly, the power law dependence 
cannot be explained only by a non-interacting percolation 
model, meaning that the statistics capture a character- 
istic property of Vanadium Oxide thin films, which may 
be a signature of self organized criticality [23( . 
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